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ABSTRACT 

We present a new technique based on multi-band near ultraviolet and optical photometry to measure 
both the stellar intrinsic properties, i.e. luminosity and effective temperature, and the interstellar dust 
extinction along the line of sight to hundreds of stars per square arcminute. The yield is twofold. On 
^Y^' , the one hand, the resulting reddening map has a very high spatial resolution, of the order of a few 

arcseconds, and can be quite effectively used in regions where the interstellar material is patchy, thus 
producing considerable differential extinction on small angular scales. On the other hand, combining 
the photometric information over a wide baseline in wavelength provides an accurate determination of 
temperature and luminosity for thousands of stars. As a test case, we present the results for the region 
\ around Supernova 1987A in the Large Magellanic Cloud imaged with the WFPC2 on board the Hubble 

Q^. Space Telescope. 

Q . Subject headings: (stars:) stars: fundamental parameters, (stars:) Hertzsprung-Russell, (ISM:) dust, 

' extinction, supernovae: individual (SN1987A), galaxies: individual (LMC) 
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1. INTRODUCTION Schlegcl et al (1998), for example, may not be satisfactory 

Qj. -1 J- r 4-1, T J i J- • j-r, TT • for detailed stellar populations studies. For example, it 

Stars provide most of the light we detect m the Um- , , i • , , ■ ; • r i i • i 

1,1 , . T, • 1 J J -1 w c cannot be used m the direction of nearby galaxies, such as 

verse and the entire chemical and dynamical evolution of ,, „ . , ,t,^oi ai ;i i ■ i i r 

„ galaxies is strongly influenced by the formation of many ^he Magellanic Clouds and M31. Also, the derived value of 

■ generations of stars. The most massive ones enrich the sur- ^"l™.'' '^^'"''^^ extremely sensitive to dust tem- 

■-' J. J. f-1-ii J 1 perature and, m any case, not accurate for high values of 

rounding medium of highly processed material and suppfy ,. ,. , a a b r-^ i 

, ^ ; r, . ^ ^ /u 4-1, --i-i, 4-1, ■ -A A extinction (Ay > 0.5 mag, Arce & Goodman 1999). More- 
vast amounts Of fcmetic energy (both With their Winds and, , 1 1 . r;, . ,, 

•4.1, 4-1, ■ 14-- 4- 1 • rp TT c over, the angular resolution of the map is pretty coarse 

even more, with their ultimate explosions as type II Su- , L i, r fr .i .ii 

X T 4- 4.-4- 4- f lu ■ ~ 6 , as compared to a lew arcseconds for the method 

pernovse). Lower mass stars constitute one of the mam ^ , i , \ , ,, .; , , i, , rr 

r -11 , 1 r-i 1 J AT-4- presented here and, thus, it does not alfow to cope etfec- 

sources of cruciaf efements, such as Carbon and Nitrogen f. , .,, , , ' ; , . , n,, • ; , i, 

/, , 14- n •JC4- ji 4- 1,1 •4-- tivefy with the patclw structure typicaf of the mterstelfar 

(through steflar winds hrst, and planetary nebufa ejection ^ t j jf 

later) in their individual evolution, and may be the main t, ■ ' , • i <■ i i ; i ; • . i ; n i . 

r 4-1 • 1 4- -r • u- i It IS certainly preferabfe to determine the amount of dust 

suppfiers of the iron group efements if, m binary systems, .„ , c ,^ . r • <• ;i 

1, J ij Tc in frout of thc stars from the fight coming from the stars 

they succeed to merge and expfode as lype la Supernovae. ,, , ttti -i . • i , , • i 

A . 1,1 4-1, 4- 4- 1 4- u- • themselves. Whife spectroscopic data may certamfy pro- 

A maior probfem that one encounters almost ubiqui- . -, , ^ . . r ,^ i i • i ;i 

, 1 1 4-j- 4- -4-1, a-4-4-i,4--4- 4-11 vidc au accurate determination of the reddening and the 

tously when studying stars is the effect that mterstelfar , -,, , . , . °„ 

dust has on the observed spectral distribution of their ^^^^If parameters, the multiplexing capabilities of current 

, , » ■ 4. r 4-1, • 4. 4. n jj ■ spectrographs are limited, at best, to a few tens of (refa- 

hght. A precise measurement of the mtersteflar reddening \ , ■ mi \i i 

, 4-1 1- r ■ 14- 4- 4-1 4- 1 4- J • 4-1, c I tivcly bright) obiects per telescope pointing, ihe method 

along the line of sight to the stars under study is the hrst / ''^ ^ ^ . . ,.° ., , . , 

, 1 • 1 , J J 4- ■ ■ 4-1, ■ 4- • ■ we describe m this paper overcomes this fimitation by us- 

step, and a cruciaf one, towards determining the intrinsic . i • • i , i • i- n 

.(• 4-1, 4- 11 14-- All 1 jj • mg onfy imaging and, thus, providing measurements of 

properties of the stellar populations. All sky reddening , , i, ; i .1 11 • 1 ;i 

, ,1 j - jf i u both the stellar parameters and the reddening along the 

maps, such as the one derived from dust temperature by ,. r • 1 , ,1 1 r , • • mi • 

fme of sight to thousands of stars per pointing. Ihe gam 

^ Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated 
by AURA, Inc., under NASA contract NAS 5-26555. 

^ On assignment from the Research and Scientific Support Department of ESA 
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in terms of statistics is obviously tremendous. As a con- 
sequence, the spatial resolution of the resulting extinc- 
tion map is extremely high, with an E{B — V) measure- 
ment every few arcseconds, even for moderately dense star 
fields. While here wc test this newly developed technique 
on HST- WFPC2 observations, the method can be applied 
to data obtained with any instrument. The quality of the 
final result, however, will obviously depend on the quality 
of the input photometry. 

The whole purpose of this paper is to describe a method 
to retrieve the stellar parameters from this kind of data 
by fitting model atmospheres to the observed fluxes. In 
section 2 we describe the dataset considered and the re- 
duction techniques adopted. In section 3 the limitations 
of the classical approach based on Color-Magnitude dia- 
grams are discussed and the method we have developed to 
overcome them via a multiband fit to the observed spectra 
is presented in section 4. The results of this procedure are 
discussed in section 5, while section 6 deals with the effects 
that the gravity and metallicity of the input model spectra 
have on the derived parameters. Finally, the spectral cov- 
erage needed to recover bot Tg// and £'(B — V) is explored 
in section 7 and the role of random and systematic errors 
is assessed in section 8. 

The broad band colors in selected WFPC2 filters and 
the corresponding extinction coefficients, which are needed 
for any application of our procedure, are provided in ap- 
pendix 9. 

2. OBSERVATIONS AND DATA REDUCTION 

Our method was developed and tested using HST- 

WFPC2 multiband observations. In particular, we have 
taken advantage of the wealth of information provided by 
the repeated observations of Supernova f 987A (SN1987A) 
and the field surrounding it in the Large Magellanic Cloud 
(LMC). Since 1994, Supernova 1987A was imaged with the 
WFPC2 every year as a part of the Siipernova INtensivc 
Study project (SINS, PI: Robert P. Kirshner). The log of 
the observations we have used is reported in Table 1. The 
WFPC2 filters do not exactly match any of the ground- 
based photometric systems and the closest analog of the 
passbands we have used is reported in the last column of 
this Table. The observations cover an almost circular por- 
tion of the sky centered on SN1987A and with a radius 
of 130", corresponding to roughly 32 pc at the distance of 
the LMC (51.4 kpc, Panagia 1999; Romaniello et al 2000). 
The true color image derived from our data is shown in 
Figure 1. 

The observations were processed through the stan- 
dard PODPS (Post Observation Data Processing System) 
pipeline for bias removal and flat fielding. In all the cases 
the available images for each filter were combined to re- 
move cosmic rays events. The long baseline in wavelength 
that they cover (more than a factor of 4, extending from 
~2,300 to ^-^9,600 A) together with the unprecedented spa- 
tial resolution of WFPC2 (~ 0".l, i.e. roughly 0.03 pc at 
the distance of the LMC) make of this dataset a unique 
and ideal one to investigate the problem of disentangling 
the effects of reddening from those of different effective 
temperatures for stars of different spectral type. 

The plate scale of the camera is 0.045 and 0.099 arc- 
sec/pixel in the PC and in the three WF chips, respec- 



tively. We performed aperture photometry following the 
prescriptions by Gilmozzi (1990) as refined by Romaniello 
(1998), i.e. measuring the flux in a circular aperture of 
2 pixels radius and the sky background value in an annu- 
lus of internal radius 3 pixels and width 2 pixels. Aperture 
photometry is perfectly adequate for our study because 
crowding is never a problem in our images. Actually, the 
average separation of stars from each other is about 1.3" 
{i.e. ~ 29 pixels in the PC chip and ~ 13 pixels in the 
WF ones) which is much larger than the WFPC2 PSF 
width. We have used the IRAF-synphot synthetic pho- 
tometry package to compute the theoretical Data Num- 
bers to be compared to the observed ones. The WFPC2 
calibration and, consequently, the synphot results are typ- 
ically accurate to within ± 5% (Whitmore 1995; Bagget 
et ai 1997). 

We identify the stars in the deepest frame (F555W, i. e. 
V band-like filter) and then measure their magnitudes in 
all of the other filters. The total number of stars identi- 
fied in this way is 21,955. If a star is not detected in a 
band, then its magnitude in that band is substituted by 
a 3a upper limit. For about 12,340 of them, i.e. 56%, 
the photometric accuracy is better than 0.1 mag in the 
F555W, F675W and F814W filters. The number of stars 
with accuracy better than 0.1 mag drops to 6,825 (31%) 
in the F439W band, and only 786 (3.6%) stars have an 
uncertainty in the F255W filter smaller than 0.2 mag. 

The brightest stars in each CCD chip are saturated (i.e. 
most stars brighter than 17.5 in the F555W passband). 
We recover their photometry either by fitting the unsat- 
urated wings of the PSF for moderately saturated stars, 
i. e. with no saturation outside the central 2 pixel radius, or 
by following the method developed by Gilliland (1994) for 
heavily saturated stars. In most cases, the achieved pho- 
tometric accuracy is better than 0.05 magnitudes. Full de- 
scription of the methods used can be found in Romaniello 
(1998). As a sanity check we have compared our photom- 
etry in the F439W and F555W bands with ground-based 
one in B and V (28 stars brighter than V ~ 20 within 30" 
from SN1987A, Walker & Suntzeff 1990) and found an 
overall excellent agreement, i.e. rms deviations less than 
0.05 mag in both filters, except for the obvious cases of 
few stars that appear as point-like in ground based im- 
ages but are resolved into more than one object on the 
WFPC2 frames. This result confirms the accuracy of our 
photometric measurements and the quality of the WFPC2 
calibration. 

We shall see in section 4 that a good photometric quality 
in different bands is needed in order to recover £^(B — V) 
and stellar temperature. In order to select stars with a 
good overall photometry we define the mean error in 5 
bands (65) {i.e. excluding the UV F255W band) as: 

_ , / "F336W ^ "F439W ^ "F555W ^ "F675W ^ "F814W 

^ (1) 

In addition to the main field described above, we have 

also used a second set of WFPC2 images pointing 8' south- 
east of SN1987A. As shown in Table 2 in this case we only 
use 4 filters and the F336W and F439W passbands are 
substituted with their broader counterparts F300W and 
F450W, respectively. Being broader, the latter are more 
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efficient and allowed us to save precious telescope time, al- 
ways a crucial issue, without losing any vital information 
(see section 7). 

3. COLOR-MAGNITUDE DIAGRAMS 

As it can be seen in Figure 1 Supernova 1987A and 
its rings arc right at the center of the field imaged with 
WFPC2 and belong to a loose cluster of blue stars (sec 
Panagia et al 2000). Also, the interstellar medium, as 
traced by the Ha emission coded in red in Figure 1, is far 
from being spatially imiform and its distribution varies on 
scales of few arcscconds (let us recall here that 1" corre- 
sponds to roughly 0.25 pc at the distance of the LMC). 
Since the spatial distribution of dust usually follows the 
one of gas (see, for example, Mathis 1990), one can expect 
the interstellar reddening to vary on similar, if not even 
smaller, scales. This is further confirmed by an inspection 
of Figure 2 where we show the Color-Magnitude diagrams 
(CMD) for different combination of bands. 

The black dots in each CMD in Figure 2 are the 6,695 
stars with ^5 < 0.1 mag, where is the mean error defined 
in equation (1). The error threshold of < 0.1 reflects 
itself as a magnitude threshold at mpsssw — 23. We esti- 
mate that the completeness at this magnitude limit is very 
close to 100% because the density of detected stars down 
to 23th magnitude is rather low, i.e. about 1 star per 4.6 
square-arcsec area (Panagia et al 2000). 

It is apparent that, despite the high quality of the mea- 
surements (internal imcertainties less than 0.1 magnitudes, 
sec also Figure 3), the various features of the CMDs, such 
as the Zero Age Main Sequence for the early type stars and 
the Red Giant clump for the more evolved populations, are 
rather "fuzzy" and not sharply defined. Although this is 
due in part to the presence of several stellar populations 
projected on each other (Panagia et al 2000), most of the 
problem arises from the fact that reddening is not quite 
uniform over the field, thus causing the points of otherwise 
identical stars to fall in appreciably different locations of 
the CMDs. 

4. THE FITTING PROCEDURE 

The large number of bands available in our dataset pro- 
vides a sort of wide-band spectroscopy which defines the 
continuum spectral distribution of each star quite well. 
Of course, the resolution is very low, of the order of 
A/AA ~ 7, but, still, important features, such as the 
Balmer jump between the HST F336W and F439W bands, 
are clearly visible. Four examples of "wide band spectra" 
are shown in Figure 4 (the flux F^ is expressed in units 
of erg cm~^ A~^) and their observed magnitudes arc 
listed in Table 3. The whole purpose of this paper is to 
describe a method to retrieve the stellar parameters from 
this kind of data by fitting model atmospheres to the ob- 
served fluxes. 

The shape of the observed spectrum is mostly deter- 
mined by the effective temperature of the star and by the 
amount of dust along the line of sight. The other quantity 
of interest, namely the solid angle subtended by the star 
(or, correspondingly, its angular radius) is obtained by a 
straight normalization of the spectral distribution. 

The general idea of our method is to recover simultane- 
ously the extinction along the line of sight and the effective 



temperature of a star by fitting the observed magnitudes 
to the ones computed using stellar atmosphere models red- 
dened by various amounts of E(B — V). The angular ra- 
dius, then, is directly obtained from the ratio of the ob- 
served flux, dcrcddencd with the best fitting E{1B — V) 
value, to the one at the surface of the star, as given by the 
best fitting model atmosphere. Thus, the two fundamental 
ingredients that go into the fit are the stellar atmosphere 
models and the reddening law: we have adopted the mod- 
els by Bessel et al (1998) and the reddening law by Scuderi 
et al (1996). 

We have chosen the models by Bessel et al (1998), as 
they provide a well tested, homogeneous set covering the 
temperature interval between 3,500 and 50,000 K in the 
wavelength range between 90 A and 160 fim. They are 
computed with an updated version of the ATLAS 9 code 
used for the Kurucz (1993) models with no allowance 
for convective overshooting, as discussed by Castelli et al 
(1997a,b). 

We have used the IRAF task synphot to compute the 
theoretical magnitudes at the surface of the star for the 
grid of temperatures provided by Bessel et al (1998) by 
convolving the model spectra with the HST+ WFPC2 re- 
sponse curves. Once the convolution is performed, we use 
a spline interpolation to get a temperature grid with a 
logarithmic step of 0.005, i.e. 1%. We have used model 
atmospheres with log(g) = 4.5 (cgs units), which is the 
typical value for Main Sequence stars, and Z= 0.3-Zq, the 
mean expected metallicity of the LMC. In section 6 we 
will analyze and discuss the effects of these two parame- 
ters on the fit and we will see that neither of them affects 
the values of the derived parameters appreciably. 

In order to faithfully model the data, we have computed 
separately the expected counts that a given spectrum pro- 
duces on each of the four WFPC2 chip detectors. Also, 
the throughput variation between decontaminations of the 
camera (Bagget et al 1996) was taken into account by spec- 
ifying the Julian date of the observation (the cont keyword 
in synphot). Details of this procedure can be found in Ro- 
maniello (1998), while model magnitudes for solar and one 
third solar metallicity are listed here in Tables (6) and (7), 
respectively. 

The other input ingredient is the reddening law. This 
is a combination of foreground extinction due to dust in 
the Milky Way and internal extinction in the LMC. We 
have taken the first component to be the same for all of 
the stars, with E{B — y)^-^ = 0.05 (see, for example, 
Schwering & Israel 1991) and we have used the Galactic 
reddening law as compiled by Savage & Mathis (1979). On 
the other hand, the amount of extinction from the LMC 
dust is one of the parameters of the; fit. We have; taken the 
LMC reddening law in the direction of our field to be the 
one determined by Scuderi et al (1996) from a study of an 
HST-FOS ultraviolet and optical spectrum of Star 2, one 
of the companions to SN1987A projected just 3".9 NW of 
it. It is, of course, crucially important to use a reddening 
law derived in the same field under study, as it may vary 
considerably from place to place on quite small angular 
scales (see, for example, Whittet 1992). These redden- 
ing laws were used to compute the extinction coefficients 
in the WFPC2 passbands for every combination of input 
model spectrum and value of E{B — V) (see Table 8). 
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The fitting procedure consists of two steps: 

• we first use a reddening-free color to determine 
whicli stars can be used to measure the reddening 
unambiguously ; 

• once this is done, we use the full information from 
the 6 bands to fit the observed spectral energy dis- 
tribution for all of the stars. 

We will now describe these two steps. 

4.1. First step: reddening-free colors 

It is a well known fact that the solution in the E{B — Y)- 
Te/ / plane may not be unique for stars later than, roughly 
speaking, AO and that different combinations of tempera- 
ture and reddening may lead to almost indistinguishable 
optical and near-UV low resolution spectra. An example 
is shown in Figure 5 where we compare Bessel et al (1998) 
models for three combinations of temperature and redden- 
ing: Teff = 5, 000 K and E{B - V) = (short-dash Uno), 
Te// = 7, 600 K reddened by E{B - V) = 0.55 (long-dash 
line) and, finally, Teff = 12,000 K and E{B - V) = 0.86 
(full line). The spectra are normalized to the F814W flux 
and crrorbars corresponding to 0.1 mag arc shown. To 
all practical purposes the spectra are identical for wave- 
lengths longer than 3,300 A and only good data in the 
ultraviolet can help to distinguish them. Unfortunately, 
though, these data are often not available, as imaging in 
the ultraviolet can only be done from space and it requires 
very long exposure times for red stars ^. The purpose of 
first step of the procedure, then, is to use a reddening-free 
color to eliminate the E{B — Y)-Tf.ff degeneracy when no 
ultraviolet observations are available. 

Let us consider the reddening- free color QubI: combina- 
tion of the magnitudes in the U, B and I passbands {e.g., 
Mihalas & Binney 1981, page 186): 



Since, by construction, the value of Qubi does not 
change with E{B — V), the effect of reddening in the 
QuBi vs (F336W - F814W) plane is to move the points 
horizontally to the right from the zero-reddening locus, as 
indicated by the arrow in Figure 6. Moreover, A(F336W — 
F814W), the horizontal distance of an observed point from 
the theoretical zero-reddening curve, is proportional to 
E(B- V): 

A (F336W - F814W) = £;(F336W - F814W) = 

- 3.3 ■ E{B - V) (4) 

whereas the observed Qubi is, to first approximation, a 
function only of the star's temperature (see Figure 6). 

As apparent from Figure 6 Qubi is not a monotonic 
function of F336W— F814W, i.e. temperature, but it is 
W-shaped'^. This, in turn, is due to the non monotonic- 
ity of the Balmer jump, the intensity of which peaks for 
AO stars and declines both for earlier and later spectral 
types (sec, for example, Allen 1973). So, stars in region 
A on Figure 6 have only one possible intersection with 
the zero reddening curve and for them the solution Teff- 
E(B — V) is unique. In the case of stars in region B, there 
are three intersections with the reddening curve, but the 
two on the right are easily rejected by noticing that, ac- 
cording to equation (4), they correspond to negative values 
of £'(B — V) and, hence, are not acceptable. The solution, 
therefore, even if not unique, is unambiguous. 

Things arc more complicated for the stars with 
F336W-F814W> 0, t.e. T,.// < 9,000 K. Here, the so- 
lution is not only non-unique, but also ambiguous. As a 
matter of fact, there are two (region C) or three (region 
D) combinations oiT^ff and E(B — V) that are absolutely 
equivalent. All of these solutions are physical and there is 
no a priori criterion the decide which one is the right one. 

The first step of the procedure is to divide the 
stars in four classes, according to their photo- 
metric error and the number of solutions in the 
Teff — E(B — V) plane. In other words, they are clas- 
sified according to the mean error in five bands ^5 intro- 
duced in equation (1) and to the number of intersections 
that they have to their left, i.e. bluer F336W-F814W 
color, with the theoretical zero-reddening curve in the 
QuBi vs. F336W— F814W plane. In order to account for 
the photometric errors, we project the measured value plus 
or minus its error. In the following, a star belongs to a cer- 
tain class if the measured point and the ±2 • a values have 
the same number of intersections. 
The classes are defined as: 

class I: stars that have only one intersection to their left. 
They are the stars in region A and B in Figure 6; 

class II: stars in region C in Figure 6. They have two in- 
tersections to their left. The solution for E{B — V) 
and Te/ / is not unique and we shall see shortly how 
we choose between the different possibilities; 

class III: stars that, in spite of their good photometry 
((55 < 0.1), have a "problematic" location in the 
QuBi vs. F336W-F814W plane. These are: 

^ For example, 4,700 seconds would be required to reach a signal to noise ratio of 5 in the F255W filter for an FOV star at the distance of the 
Large Magellanic Cloud. 

3 The additional hook at F336W-F814W> 3.5 is caused by the WFPC22 F336W filter red leak. 



OUBI = (U - B) 



E{\J - B) 



E{B - 1) 

QuBi is reddening-free in the sense that: 



(B-I) 



(2) 



Q 



UBI 



(U-B)„ + i?(U-B) 



E(V-B) 



(3) 



E{B - I) 

[(B-I)o + £;(B-I)] -QuBLo 
where the subscript "0" indicates the unreddened quan- 
tities. QuBi is not exactly equal to Qubi.o because the 
coefficient in equation (2) is not a constant, but depends 
on the temperature of the star and on the amount of red- 
dening in front of it (see Appendix 9). With the redden- 
ing law of Scuderi et al (1996), and using the WFPC2 
filters F336W (U), F439W (B) and F814W (I), one 
finds i;(F336W-F439W)/£;(F439W-F814W) ~ 0.42 
for Teff > 5, 500 and this is the value will use. 

The reddening-free color Qubi computed with the mod- 
els by Bessel et al (1998) for the WFPC2 filters is plotted 
in Figure 6 versus (F336W — F814W). The curve is for 
E{B — V) = and the arrow is the reddening vector for 
E{B - V) = 0.2 mag. 
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• stars whose horizontal projection never inter- 
sects the zero-reddening curve to their left; 

• stars that belong to region D Figure 6 and, 
hence, have three intersections to their left. 
As we shall discuss in section 6, in this re- 
gion the particular combination of colors we 
have chosen is sensitive to metallicity and this 
classification scheme is not reliable; 

• stars that fall above the peak at 
F336W-F814W= 1, roughly along the long- 
dash line in Figure 6, for which the F336W 
flux is too strong as compared to the one in 
the other bands. As discussed in Panagia et 
al (2000), they are UV excess objects, most 
likely T Tauri stars and the U-band excess is 
due to the presence of a circumstcUar disk ac- 
creting onto the star. As a consequence, their 
spectrum is not well modeled by the one of a 
normal photosphere and no reliable iS(B — V) 
can be derived by comparing their colors to 
those of a normal photosphere. 

class IV: finally, stars with poor overall photometry, i.e. 
^5 > 0.1, are assigned to this class, independent of 
the number of intersections. 

As an example, the result of the division into classes is 
shown in Figure 7 for the WF3 chip of the July 1997 ob- 
servations. Class I stars are hotter than about 10,000 K, 
whereas class II stars arc between 6,750 and 8,500 K. In 
the particular case we are discussing here, out of a total 
number of 3,139 stars, 40 belong to class I (1.3%), 388 to 
class II (12.4%) and 615 to class III (19.6%). The rest of 
them (2096 or 66.8%) belong to class IV. 

Again, for every star the measured Qubi gives a first 
order estimate of the temperature, see Figure 6, and the 
excess in £^(F336W — F814W) gives a first order estimate 
of i?(B — V) from equation (4). These are, then, used as 
an input in the second step of the procedure, the 6-band 
fit, which we will now describe. 

4.2. Second step: 6-band fit 

In the first step we have used only a subset of the 6 
bands available. Now that the number of possible solu- 
tions in the Te//--E(B — V) plane is assessed for every star, 
we can proceed to fit the entire spectrum at once with a 
technique using the information gathered in the first 
step as a starting point for the fit: 

class I: The fit is performed in a neighborhood of 
±0.1 mag around the E{B — V) value deducted 
from the projection in the Qubi vs. F336W-F814W 
plane; 

class II: We choose the smallest of the two E{B — V) 
values from the projection described above as the 
starting point and perform the fit in a neighborhood 
around it, defined as to avoid the other possible so- 
lutions. This is somewhat arbitrary, but an inspec- 
tion of Figure 7, for example, does not indicate any 
compelling evidence for a large number of stars with 
extinction much larger than the mean value. Let us 
stress here that the reddening is still computed in- 
dividually for each star; 



class III and IV: for the various reasons mentioned 
above, photometry is not reliable enough to allow 
to solve for all of the three unknowns simultane- 
ously and the reddening is derived from class I and 
II stars. Wc have explored two options: either using 
the mean value of the 4 closest class I and II neigh- 
bors or the mean over the entire field. In the case of 
the field of SN1987A the result is the same, as the 
extinction does not show any significant spatial pat- 
tern and the local mean is mostly equal to the gen- 
eral one. However, when the reddening shows ap- 
preciable small scale correlations, however, the local 
mean does give better results, most notably a nar- 
rower Main Sequence, and the first option should be 
preferred (see, for example, the case of NGC 6822 
in Bianchi et al 2001). 



The details of the fit a^^e described in Appendix 9. 
While in the selection process we have used constant ex- 
tinction coefficients to build Qxjbi (sec equation (2)), in the 
6-band fit they vary according to the Tg// and i?(B — V) 
values (see Table 8 in Appendix 9). 

This procedure assures total control over the fitting pro- 
cess. In fact, even in the case of multiple intersections, it 
provides a criterion to choose the most reasonable solution. 
Once this is done, the 6 band fit uses all of the informa- 
tion available and, thus, gives the best possible answer. 
Moreover, the whole procedure is very eSicient in terms of 
computer CPU time usage, since it limits the size of the 
grid in the 6 band fit, which is the most time-consuming 
part of it all, only to a small region in the parameter space 
around the physical minimum in the Teff-E{B — V) plane. 



5. RESULTS OF THE FIT 

As an example of the results of the procedure, in Fig- 
ure 8 we show the spectra of the same four stars of Figure 4 
together with the best fitting models. The derived param- 
eters and their errors are listed in Table 4. 

The HR diagram for the 21,955 stars in the field of 
SN1987A is displayed in Figure 9. It is interesting to com- 
pare it to what one would have obtained from a "blind" fit 
of the spectrum, i.e. by skipping the selection described in 
section 4.1 and performing directly the 6-band fit of sec- 
tion 4.2. This is shown in the left panel of Figure 10. 

The fit for the stars marked in grey in the left panel of 
Figure 10 is obviously wrong as they occupy an impossi- 
ble region of the HR diagram. An inspection of the right 
panel of the same Figure reveals that all of these stars, and 
only them, were fitted with too high a value of E(B — V) 
and, consequently, of Teff- These are all stars that fall 
in regions C or D in Figure 6 for which the broad band 
spectrum can be interpreted either as a heavily reddened 
hot star or as a cold one with little dust in front of it. The 
fit just selects the solution with the lowest x^) even if it 
is marginally lower than other possible minima. Which 
one of the possible solution has the lowest depends on 
such things as the E{B — V) and Te// grids used in the 
fit and, in the case shown in Figure 10, the wrong one is 
chosen for two stars out of three. For many of these stars 
the available photometry is limited to F439W and redder 
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filters. As we will see in section 7, then, the data are very 
little sensitive to temperature and these cold stars can be 
mistaken for objects as hot as 50,000 K. 

On the other hand, selecting the right solution before 
performing the fit eliminates the problem almost com- 
pletely. There still are roughly 600 stars (less than 3% 
of the total) under the Zero Age Main Sequence locus. 
However, they are mostly class IV objects, i.e. stars with 
very poor photometry, or stars with peculiar colors, i.e. the 
T Tauri stars described in detail by Panagia et al (2000) 
or unresolved binaries. 

The stars with multiple solutions are the late-type ones. 
Being able to measure directly the reddening for them is 
important for two reasons. Firstly, as low mass stars are 
much more numerous than higher mass ones, the num- 
ber of reddening determinations increases dramatically. In 
fact the late type stars for which eventually we are able 
to measure the reddening are typically a factor of 10 more 
numerous than those earlier that AO. Secondly, stars of 
different masses have different lifetimes and there is no a 
priori reason why different generations of stars should be 
affected by the same amount of extinction (see Zaritsky 
1999, for an example in the LMC itself). By extending 
our analysis to a wide range of masses we can check for 
population-dependent effects and deredden the different 
generations of stars in the most appropriate way. 

In addition to measuring the intrinsic properties of the 
stars shown in Figure 9, our procedure also provides indi- 
vidual E{B — V) measurements for more than 2,500 stars, 
on average one every 13 square arcseconds. The resulting 
histogram and spatial distribution of interstellar redden- 
ing are shown in Figure 11 (see also Panagia et al 2000). 
As it can be seen, the initial suspicion is c:oiifirnicd and 
the reddening is indeed very patchy and varies on scales 
of a few arcseconds. 

Besides providing a very accurate reddening map, mea- 
suring individual values of E{B — V) also leads to a nar- 
rowing of the various features in the Color-Magnitude di- 
agram, in particular the Main Sequence. This is of great 
help when interpreting the data, because it allows one to 
disentangle the broadening due to differential extinction 
from the one caused, for example, by the overlap of sev- 
eral generations of stars or by the presence of binary stars. 

In Figure 12 we compare the stars in the 
F439W vs (F439W-F814) diagram before and after red- 
dening correction. The density contours in panels (a) 
and (b) are spaced by factors of 2. Apart from the ob- 
vious fact that the corrected stars are bluer and more 
luminous, it is apparent that the features in the reddening 
corrected diagram are narrower than the corresponding 
ones in the observed CMD. In order to quantify this effect 
let us consider a color cut through the Main Sequence, 
like the one displayed in panel (c) of Figure 12 (see also 
Figure 13). Here, the blue line histogram represents the 
color distribution of the dcreddened Main Sequence stars 
with ^5 < 0.1 in the 21 >m(F439W)o > 22 magnitude 
interval. The red line histogram gives the distribution of 
the same stars before reddening correction. To facilitate 
the comparison, this latter is shifted in color so as to have 
the same mode as the dereddened one. As apparent, the 
original distribution is definitely more spread out in color: 
its rms is 24% larger than the corrected one (0.21 mag 



versus 0.17). This is the case all along the Main Sequence, 
as shown in Figure 13, where we plot the color distribu- 
tion for all the stars with 5 5 < 0.1 in different magnitude 
bins. The dereddened histograms (blue line) always dis- 
play a smaller scatter than the observed ones (red line), 
confirming the effectiveness of the correction procedure. 

As it can be expected, the color distribution is even nar- 
rower when only the stars with individual reddening de- 
terminations, i.e. those we have classified as class I and II, 
are considered. This is shown in Figure 14. In this case the 
rms of the reddening-corrected Main Sequence histogram, 
again in the 21 >m(F439W)o > 22 magnitude interval, is 

0. 14 mag, while the one of the same stars before correc- 
tion is 0.18 mag, i.e. 29% larger. A simple Kolmogorov- 
Smirnov test shows that the widths of the two distributions 
are different at the 99.9995% level. 

To conclude, let us remark again that Class III and IV 
stars were not corrected with and individually determined 
value of i?(B — V), but, rather, with the mean value from 
2,510 Class I and II neighbors. It is important to realize 
that, given the large number of these stars, the mean red- 
dening used to correct Class III and IV stars is on average 
very accurate. Of course, it is still possible that a few class 
HI and IV stars have, in reality, extinction values signifi- 
cantly different from their local mean, but the global effect 
is negligible. 

The full discussion of the stellar content of the region 
around SN1987A can be found in Romaniello (1998) and 
Panagia et al (2000). 

6. THE ROLE OF GRAVITY AND METALLICITY 

So far we have used model atmospheres only with one 
surface gravity (log(.g) = 4.5, appropriate for Main Se- 
quence stars) and metallicity (Z= 0.3 • Zq, the expected 
mean value for the LMC). Even though these are sensible 
assumptions for the environment under study, we have to 
verify their influence on the results we derive. 

Let us start by considering the surface grav- 
ity. The theoretical zero-reddening locus in the 
QuBi vs. F336W— F814W plane is shown in Figure 15 for 
three values of surface gravity, namely log{g) — 4 (dotted 
Hue), 4.5 (solid line) and 5 (dashed line). As it can be seen, 
the three curves overlap perfectly for F336W— F814W< 
-0.5, i.e. Teff > 11,000 K and F336W-F814W> 1.5, 

1. e. Teff ^ 5,500 K. Elsewhere, however, different surface 
gravities lead to different expected zero-reddening loci in 
the QuBi vs. F336W-F814W diagram. In particular, this 
influences the division in classes that is the foundation of 
the method described above. 

However, as it can be seen, for example, in Figure 14, 
only stars on or near the Main Sequence are used to de- 
termine the reddening. Adopting the stellar evolutionary 
models by Brocato & Castellani (1993) and Cassisi et al 
(1994), we see that for Z= 0.3 • Zq the gravity is nearly 
constant along the Main Sequence and that its value is 
log (5) =4.5, the value we have used in the fit. This is 
a sound assumption over a large range of metallicities be- 
cause extensive model calculations show that surface grav- 
ity variations as a function of chemical composition are 
quite smaU. For example, at (F336W - F814W)o = 
the surface gravity ranges from log{g) — 4.4 for a solar 
metallicity star to log{g) ~ 4.6 for Z=Z0/2O. Considering 
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that for every metallicity the gravity is almost constant 
along the Main Sequence one can safely adopt the value 
of log(g) = 4.5, independent of chemical composition. 

Once a star evolves off the Main Sequence its surface 
gravity changes appreciably. For example, \og(g) ~ 3 is 
a typical value for the stars in the Red Giant clump, i.e. 
the stars at log(Te|y ) ~ 3.7 and log(L/L0) ~ 1.75 in Fig- 
ure 9. Even though the Qubi vs. F336W-F814W relation 
is, indeed, sensitive to gravity, the result of our global fit 
is not. This is because, as we have seen, we do not fit si- 
multaneously temperature and reddening for stars off the 
Main Sequence, but rather we assign to them the mean 
E{B — V) value of their class I and II neighbors, which are 
on the Main Sequence and, hence, have log(.g) ~ 4.5. For 
stars that are off the Main Sequence, we just perform the 
fit to solve for the radius and effective temperature. Once 
the reddening is correctly determined, the result of the 
multi-band fit is no longer sensitive to the surface grav- 
ity adopted in computing the theoretical colors. We have 
checked this by dividing the stars into classes using the 
QuBi vs. F336W - F814W relation for log{g) = 4.5, and 
then fitting the stars with colors computed for various val- 
ues of the gravity. The maximum variations in tempera- 
ture for the stars on the Red Giant branch are of 3%, i.e. 
much smaller than the error deriving from the fit itself, 
and with no systematic trend with gravity. Given this re- 
sult, we have decided to use \og{g) = 4.5 all throughout 
our analysis. 

In addition to influencing the stellar structure, the 
chemical composition also plays a role in determining the 
star's spectrum, through the atmospheric opacity. As a 
consequence, the colors also depend on this parameter. 
The OuBi vs. F336W - F814W relation for different val- 
ues of Z is shown in Figure 16. Again, as it can be seen, 
the variations are negligible for the stars for which the red- 
dening is an individual fit parameter, i.e. Class I and II 
stars. 

Also in this case, then, once the subdivision into classes 
is made and the local mean reddening for Class III and 
IV stars is determined from the individual values of Class 
I and II neighbors, the result of the multi-band fit is in- 
sensitive to the details of the stellar atmospheres. The 
comparison of the fits performed with colors computed for 
different metallicities shows that the derived parameters 
do not depend appreciably on metallicity. 

We conclude that in applying our fitting procedure, as 
described in section 4, one can safely use colors computed 
with models by Bessel et al (1998) for Z= 0.3 • Zq and 
log{g) = 4.5 to derive accurate stellar parameters and in- 
terstellar absorptions for all stars. 

7. THE BANDS NEEDED 

For every star there are three unknown quantities: the 
angular radius, the effective temperature and the redden- 
ing. Therefore, for each star one needs at least data in 
three bands to make the problem mathematically well de- 
fined. For good results the measurements in those three 
bands must also be rather accurate, thus requiring that 
the bands be selected so as to fall around the wavelength 
of the maximum of the spectral distribution of any given 
star, and have a suitable wavelength spacing to assure a 

* The hook at very low temperatures is due to the red leak in the V 



sufficient baseline for the spectral shape analysis. In prac- 
tice this is achieved with bands separated from each other 
by about 30-40% of their central wavelengths. Moreover, 
when dealing with a stellar population one wants to study 
stars within a wide range of temperatures, and, therefore, 
more than three bands are absolutely needed. In this sec- 
tion we will investigate how many and which bands are 
needed to constrain the fit. 

Let us start by noticing that a U-like band is mandatory 
to recover both Teff and i?(B — V). To illustrate this point 
let us define the reddening-free color Qbvi combination of 
B, V and I. By analogy with equation (2): 

gBvi-(B-V)-^^.(V-I) (5) 

Qbvi is shown in Figure 17 as a function of F439W — 
F814W in panel (a) and Teff in panel (b). 

Since the the B, V and I filters do not include the Balmer 
jump, the relation between Qbvi and B— / is monotonic"', 
thus solving the problem of ambiguities in the solution. 
However, one immediately notices that the slope of the 
curve is much shallower than in the Qubi case. In fact, 
the total excursion of Qbvi between 3,500 and 50,000 K 
is only 0.4 mag and it is almost totally insensitive to tem- 
peratures higher than 6,000 K, i.e. F439W - F814W < 1. 

The combined effect of the weak dependence of Qbvi 
on temperature and of photometric errors is to make ev- 
ery combination _E(B — V)-T(.// along the reddening vec- 
tor essentially equivalent. This is shown in Figure 18 
and Figure 19 where maps of the fit using only four 
bands (F439W, F555W, F675W, and F814W, left panel) 
are compared to those using all six of them (right panel) . 
As it can be seen, when only 4 bands are used the has 
no distinct minima, and, hence, it is not possible to deter- 
mine both the reddening and the temperature neither for 
hot stars for which the solution is unambiguous (Figure 18) 
nor for the cold ones with multiple minima (Figure 19). 

Also, let us notice here that, while observations at wave- 
lengths shorter than the Balmer jump are mandatory to 
recover both Tgff and E{B — Y), it is preferable not to 
go too far in the ultraviolet, unless one can afford to make 
observations in more than one UV band. Qualitatively, 
the reason for this is that shorter wavelengths are more 
affected by interstellar reddening and small differences in 
E(B — V) along different lines of sight may result in large 
differences in the observed flux, thus making the ultravi- 
olet flux of a moderately reddened hot star almost indis- 
tinguishable from the one of an unreddened, colder one. 
Quantitatively, the situation is illustrated in Figure 20, 
where we plot the la contour levels for artificial stars of 
different effective temperature as observed through 4 opti- 
cal filters (F439W, F555W, F675W and F814W) and one 
UV one: F170W (dot-dashed line), F255W (dashed line) 
or F336W (full line). The contours, computed assuming 
the same photometric error in all 3 UV filters, overlap 
almost perfectly. This means that one does not gain in ac- 
curacy by observing at shorter wavelengths, even for hot 
stars. Since the sensitivity of normal CCDs and the emis- 
sion of cold stars both drop dramatically the further one 
moves to the ultraviolet, a filter in the region of the U 

'2 F439W filter. 
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band is the ideal choice for the filter bluer of the Balmer 
jump. 

As we have pointed out at the beginning of this section, 
there are 3 unknowns in the fit so that 6 bands are actually 
redundant and 4 could be enough to cover the wide range 
of temperatures spanned by stars in regions with mixed 
stellar populations. In fact, the method works rather well 
also with only 4 bands, provided that they encompass 
the Balmer jump. This is shown in Figure 21, where the 
HR diagram for the control field located 8' south-cast of 
SN1987A, derived following our procedure, is shown. In 
this CclSG^ clS reported in Table 2, the four available photo- 
metric bands were F300W, F450W, F675W and F814W. 
The star symbol shows the location of the brightest star 
in the WFPC2 field of view. It is so heavily saturated 
in the WFPC2 chip that its image was spilling over sev- 
eral neighboring stars thus making impossible to make di- 
rect measurements, and, therefore, its photometry is taken 
from Fitzpatrick (1998). 

As it can be seen, the observed Main Sequence, although 
still in very good agreement with the theoretical one, is 
less sharply defined than the one in the field of SN1987A 
(see Figure 9) and more stars fall under the theoretical 
ZAMS. Using a reduced set of filters has affected the re- 
sult in two ways. Firstly, there is only one band blueward 
of the Balmer jump (F300W) instead of two (F255W and 
F336W): the fit is less sensitive to temperature for hot 
stars and, hence, for them, the results are less accurate. 
Secondly, since the observations were taken as "parallels" 
to FOS primaries, we did not have enough time to make 
observations both in the F555W, or a similar band like 
F606W, and in the F675W filter. We opted for this latter 
one to identify stars with Ha emission more accurately. As 
a result, though, the photometric points are not spaced in 
wavelength quite evenly. As a consequence, the fit for the 
cooler stars, the ones that are not well exposed in F300W 
and F450W bands, has to rely on a short wavelength base- 
line and, thus, is less precise than it would have been if 
observations at roughly 5000 A had been available. This 
also results in a larger fraction of stars below the Main 
Sequence (roughly 600 out of a total of nearly 13,100, i.e. 
4.6%) when compared to the main field (roughly 600 out 
of almost 22,000, i.e. less than 3%). In general, houjever. 
the resulting stellar parameters are still of very good qual- 
ity, with more than 30% of the stars with an error in Tg/ / 
smaller than 10%. 

8. THE EFFECT OF RANDOM AND SYSTEMATIC 
PHOTOMETRIC ERRORS 

In order to assess the intrinsic precision of the procedure 
described in this paper we have tested it on simulations. 
Artificial stars were created with E(Q — V) = 0.2 in the 6 
filters listed in Table 1 , evenly distributed in the logarithm 
of effective temperature between 3,500 and 50,000 K, the 
range covered by Bessel et al (1998) models. They were, 
then, fitted according to the prescriptions of section 4. 

8.1. Random errors 

The effects of random errors were mimicked by creat- 
ing three groups of 10,000 stars and adding to the model 
magnitudes errors drawn from Gaussians of different rms 
(0.02, 0.05 and 0.1 mag). These stars were then fitted 



using only 4, 5 or the complete set of 6 bands. 

The situation for a five-band fit (F336W, F439W, 
F555W, F675W and F814W) and varying the photomet- 
ric precision is shown in Figure 22. There we plot the 
ratio of the output to input temperature as a function 
of the input temperature for different values of the error 
rms: 0.1 mag in panel (a), 0.05 mag in panel (b) and 
0.02 mag in panel (c). In each panel the photometric er- 
rors were drawn from the same Gaussian, independent of 
the star's temperature. As it can be seen, there is no sys- 
tematic trend and the points are mostly scattered around 
Tout = Tin- The values of the mean and the dispersion 
are reported in Table 5. 

The biggest deviation from the random scatter is rep- 
resented by the two features at log(Tin) — 3.95, which 
are most prominent in panel (c). These are stars close to 
the turning point at F336W - F814W ~ in Figure 6 
and were assigned the wrong value of E(Q — V) during 
the selection process which constitutes the first step of the 
procedure. Note that this effect depends on the input pho- 
tometric error because the class assignment is based on it. 
In the worst case, when the error is 0.02 mag, however, this 
systcmatics affects less than 10% of the stars at that tem- 
perature, split almost evenly between stars with too high 
and too low output temperature. This effect is negligible 
for studies of stellar populations. However if one is deal- 
ing with special categories of objects, for example variable 
or otherwise peculiar stars, which represent a small frac- 
tion of the total population in the 3.92 < log(Te//) < 4 
range one should be aware that random errors may create 
false outliers with temperatures deviating as much as 20% 
for up 10% of the stars in this temperature range. For 
this type of studies it is safer not to determine extinction 
for individual stars but, rather, to adopt the value of the 
neighbors not affected by this problem. 

It is interesting to notice in all of the panels in Fig- 
ure 22 that the dispersion increases for hotter stars. This 
is because broad band colors become progressively more 
insensitive to temperature as it increases. In fact, the 
scatter at high temperatures becomes smaller when the 
F255W filter is added to the fit, as shown in Figure 23. 
Of course, in reality hotter Main Sequence stars are also 
brighter than colder ones and, thus, have smaller photo- 
metric error which compensate for the lower sensitivity of 
colors to high effective temperatures. Again, see Table 5 
for the values of the mean and dispersion. 

8.2. Systematic errors 

In addition to random errors that vary from star to 
star and from band to band the measured fluxes are also 
affected by systematic errors in the calibrations of each 
band, i.e. errors in the photometric zcropoints. These, of 
course, are different from band to band, but the same for 
all the stars in each band. Again, we have explored three 
cases (0.02, 0.05 and 0.1 mag) and the results are shown 
in Figure 24 that displays the distributions of the tem- 
perature deviations log(Tout/Tin) for four representative 
temperatures, i.e. 5,000, 10,000, 20,000 and 40,000 K. 

From an inspection of Figure 24 one notices again that, 
given the same error, the temperature measured for hot 
stars is intrinsically more uncertain than for cold ones, 
with the dispersion growing from less than 1% for stars at 



Accurate stellar population studies 



9 



5000 K, to as much as 12% at 40,000 K. The sharp peak at 
log(Tout/Ti„) ~ 0.1 found for Tin ^ 40, 000 K (Figure 24, 
panel (d)) is populated by stars for which the fit, starting 
from Tin = 40, 000 K, would indicate a temperature higher 
than 50,000 K, the highest in the Bessel ct al (1998) grid. 
Even though the peak is just an artifact, it is clear that 
calibration errors much smaller than 0.1 mag are necessary 
for the data to be useful at all temperatures. 

9. SUMMARY AND CONCLUSIONS 

In this paper we have presented and discussed a new 
technique that uses multiband photometry to measure 
reddening, effective temperature and, given the distance, 
bolometric luminosity for resolved stars. The method is 
based on broad band optical and near ultraviolet pho- 
tometry and, in the test case presented here for the field 
around SN1987A in the Large Magellanic Cloud, yields 
to a reddening determination every 13 square arcseconds: 
a fine grid indeed! Being able to measure the reddening 
to thousands of individual stars is of great importance in 
all the astrophysical environments in which the interstel- 
lar medium is clumpy and, therefore, produces differential 
reddening, as is the case, for example, of virtually every 
star forming region. 

There are 3 unknown quantities in the fit: the radius, 
the effective temperature and the reddening. Formally, 
then, at least three photometric points are required for 
each star to make the problem well defined, and more are 
needed when dealing with stellar population encompass- 
ing a wide range of temperatures. Obviously, the more 
bands are available, the better the result of the fit. As we 
have discussed thoroughly, at least one of the bands has to 
fall at wavelengths shorter than the Balmer jump because 
the optical fluxes alone are not enough to reliably recover 
both Te// and £(B — V). On the other hand, though, 
in the case in which observing time constraints limit to 
one the the number of bands that can be used for UV 
observations, it is not necessary nor advisable to resort 
to bands at extremely short wavelengths, and observing 
just blueward of the Balmer jump, a U-like band, is the 
best choice. Of course, observing in more than one ultra- 
violet band can only help in reducing the errors on the 
derived quantities. Although we have used the WFPC2 
counterparts to the classical Johnson-Cousin passbands to 
illustrate our method, it proved to work almost equally 
well when broader filters were used. Broader filters means 
more photons detected per unit time and, ultimately, a 
more efficient use of the telescope time at one's disposal. 

Once the stars are accurately placed in the Hertzprung- 
Russcl diagram, the physical properties of the stellar pop- 
ulation, such as star formation history, mass function etc, 
can be studied in detail. In this respect, let us stress 
again the importance of determining the reddening from 
the same stars which are the ultimate target of the investi- 
gation, without having to resort to uncontrollable assump- 
tions on the magnitude and distribution of the interstellar 
extinction. 

To conclude, we refer the reader to Panagia et al (2000) 
for a discussion of the properties of the young stellar pop- 
ulation in the field around SN1987A and to Romaniello 
et al (2000) for an application of the dereddend Color- 
Magnitude diagram to measuring the distance to the Large 



Magellanic Cloud. More work on the derivation of the Ini- 
tial Mass Function and the Star Formation Rate in the SN 
1987A field is in progress, and the results will be presented 
in a forthcoming paper Romaniello et al (2002) . 
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APPENDIX A: THE FITTING PROCEDURE 

As repeatedly stated, the general idea of our method 
is to recover simultaneously the extinction along the line 
of sight and the effective temperature of a star by fitting 
the observed magnitudes to the ones computed using stel- 
lar atmosphere models reddened by various amounts of 
E{B — V). In this Appendix we describe the details of the 
fitting procedure. 

Among the three parameters we want to determine, only 
two change the shape of the spectrum, i.e. E{B — V) and 
effective temperature (Te//). The third one, i.e. the ra- 
dius i?, causes just a rigid shift of the fiux on a logarithmic 
scale. So, we first determine the shape of the spectrum by 
solving simultaneously for E{B — V) and Te//, then we 
compute the radius of the star. For every star, the follow- 
ing steps are performed: 

1. The input spectrum is normalized to the band 

with the highest photometric accuracy, usually 
the F555W, to eliminate the overall normalization, 
which depends an the star's radius. Obviously, the 

theoretical spectra are normalized in the same way. 
Let mobs^i be the normalized observed magnitude in 
filter i and ffi",, the normalized model magnitude 
for temperature j — th in the same filter; 

2. The theoretical magnitudes are reddened, for var- 
ious values of E{B — V). If TZi is the extinction 
coefficient in the band i due to the LMC dust and 
'^i,MW is the one due to our Galaxy: 



m. 



y [E{B - V)] = m^,,,,. + [7^. • EiB - V)] 



[ni,MwE{B-V)^^] (6) 
with £^(B — V)]y[.y^ = 0.05, as already discussed. 
The extinction coefficients for selected WFPC2 fil- 
ters are listed in Table 8; 

3. For any given temperature and reddening, we have 
computed the as: 



^ ^ ''rrims^ij [EiB - V)] - mobs,i'^'^ 



i=l 



(7) 
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where N is the number of photometric points and 
Smobs,i is the error on the observed magnitude 
i^obs,i- The quantity 5mobs,i^ of course, is the sum 
in quadrature of the photometric error on the mag- 
nitude in band i phis the one on the magnitude in 
the band to which the spectrum was normaUzed. 
The index i in the sum runs on all bands, but the 
one that was used to normalize the spectrum. If a 
star is not detected in a certain band, typically in 
the UV, a 3(t upper limit for the flux is used. The 
models that predict a flux in that band higher than 
this limit are rejected. 

Of course, the best-fit model corresponds to the 
combination temperature, and E{Q — V) that min- 
imizes x| [^(B - V)] ; 

4. Once the minimum x| is found, the errors in temper- 
ature ((5Te//) and reddening {5E{& — V)) are com- 
puted from the map; 

5. Now that the shape of the spectrum is determined, 
we compute the radius of the star as the error- 
weighted mean of the logarithmic shifts required to 
match the model to the observed flux in every band. 
Since the theoretical fluxes are flux densities at the 
surface of the star, the angular radius p derived from 
band i is: 



log(pi) = -0.2- (ml_,-m^^*_,) 



(8) 



where m{^* ^ is the magnitude in filter i of the best- 
fitting theoretical spectrum and mjjj,^ ^ is the ob- 
served magnitude, dereddened with the the best- 
fitting value of £'(B — V). The best estimate of the 
angular radius p, then, is: 



En I ,-o 



(9) 



where n is the number of bands with actual flux 
measurements. The associated error 5p, then, is: 



may introduce additional scattering in the HR dia- 
gram. 

For example, in the case of the LMC, we have 
adopted the value originally measured purely ge- 
ometrically by Panagia et al (1991) and later im- 
proved by Panagia (1999). Using the rings around 
SN1987A, they found D = 51.4 ± 1.2 kpc; 

6. Finally, the intrinsic luminosity (L) can be com- 
puted according to the black body law: 



L 



R 
R0 



(13) 



where Te//,© is the temperature of the Sun. 

Since the errors in temperature and radius are cor- 
related, the error on the luminosity cannot be com- 
puted in a straightforward way from the previous 
equation. It is more convenient to write the lumi- 
nosity in another, completely equivalent form. Let 
us flrst define the quantities: 



I, 



■bs.: 



i.e. the observed flux in filter i, and: 



b = 



fit 



(14) 



(15) 



where ^ is the flux in filter i of the best-fitting 
model. The quantity b is the equivalent of the bolo- 
metric correction, relating the total flux in n bands 
to the bolometric luminosity per unit surface crT^j^. 
We can now write the luminosity as: 

L = AttD^ lobs h (16) 

where lobs = J2i hbs,i- AH the quantities in this 
last equation are independent from each other and, 
hence, it is straightforward to compute the error on 
the luminosity: 



{5p) 



(n - 1) • X;r=i l/S'^mobs,i 



(10) 



To convert the angular radius of a star to the linear 
one (i?) , it is necessary to know its distance D: 




+ 



bs 



'6h 



(17) 



R = p- 



D 

Rn 



(11) 



in this case, R is expressed in units of the solar ra- 
dius R0. 

Given the uncertainty on the distance, the error on 
the linear radius {5R) is: 



(12) 



This error does not take into account any intrinsic 
dispersion in the distance of individual stars, which 




Slobs 



i^obs 



and: 



\ i=l 



'-l)-h,obs] (18) 



Sb 

T 



db 



dT. 



6T. 



eff 



eff 



(19) 



In this last equation, it is clear that the first term 
on the right-hand side, db/dT^ff, is completely de- 
termined by the models, whereas the second one, 
STeff, is the error on the temperature from the fit 
computed in step 4. 
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In the case of the LMC, 5D/D = 0.02 (Panagia 
1999). 

Let us stress once again that the stars may not be 
all at the same distance and that this may intro- 
duce additional scattering in luminosity. From pho- 
tometry alone, distances can be measured only in a 
statistical sense for an entire population and not for 
individual stars. However, since the intrinsic depth 
the LMC is estimated to be smaller than 600 pc 
(Crotts et al 1995), this uncertainty is negligible. 

APPENDIXB: BROAD BAND COLORS AND EXTINCTION 
COEFFICIENTS IN SELECTED WFPC2 FILTERS 

The WFPC2 is equipped with a large set of broad, 
medium and narrow band filters for imaging (see Biretta 
et al 2000, for a complete list and full details). Some of 
them closely match filters in the standard Johnson-Cousin 
(e.g. the F555W and the V, the F439W and the B, etc) 
or Stromgren (e.g. the F41GM and the v, the F467M and 
the b, etc) systems, while some other have no direct coun- 
terpart in any standard ground-based set, e.g., the F450W 
and F6G6W, just to name two, or, obviously, all ultraviolet 
passbands. 

Color transformations have been devised to convert the 
WFPC2 magnitudes into the corresponding ones in the 
Johnson-Cousin system, but, they are limited to only a 
handful of passbands. Moreover, quoting directly from 
Biretta et al (2000), they "should be used with caution for 
quantitative work" , because they are just first order ap- 
proximations and unavoidably degrade the quality of the 
original photometry. A much better option, when compar- 
ing WFPC2 photometry to theoretical expectations, is to 
compute these latter ones directly in the passbands used in 
the observations. As explained in section 4, this is what we 
have done in this paper, thus minimizing the uncertainties 
in the comparison between the observed WFPC2 magni- 
tudes and the synthetic ones computed using the Bessel et 
al (1998) models. 

In Tables 6 and 7 we list the magnitudes in widely used 
WFPC2 passbands for solar metallicity and for [Fe/H] = 
—0.5, the appropriate value for the LMC, respectively. 
They were computed using the IRAF-synphot synthetic 
photometry package and, for the sake of generality, the 
magnitudes in the ultraviolet filters are referred to the de- 
contaminated camera. The throughput degradation due 
to the accumulation of material between decontamination 
can be easily computed with the formula; given in Bagget 
et al (1996). The magnitudes are for a 1 R© star as seen 
from a distance of 10 pc and scaling them for any radius 
and distance using equation (20) is straightforward. 

The other important quantities that one needs in or- 
der to interpret the data are the extinction coefficients. 
Again, it is very important to compute them in the very 
same passbands used for the observations. The values of 
the normalized extinction TZx = A\/E{B — V) for the fil- 
ters of Tables 6 and 7 computed for the galactic reddening 
law as compiled by Savage & Mathis (1979) and the models 
atmospheres by Bessel et al (1998) are listed in Table 8 for 
various values of the effective temperature and the color 
excess i?(B - V). 

Using the entries of Tables 6, 7 and 8, then, it is very 
easy to model a star of metallicity Z, surface gravity g, 



effective temperature Tg/j, radius R (in solar units) as 
seen at a distance D (in parsecs) through a layer of dust 
characterized by a reddening i5(B — V): 



m(Z, g, Teff, R, D, E{B - V)) = mo(Z, g, T^fj) + 



-5 



-I- 5 • log 



+7^(Te//,i?(B-V))•i^(B-V) (20) 

where mo(Z, (/, Tg//, .^) is the appropriate value, ei- 
ther read or interpolated from Table 6 or 7, and 
7^(Te//, £;(B - V)) can be taken from Table 8. 

Ultraviolet WFPC2 filters are notoriously affected by 
red leak, and a quitc^ substantial one when red stars are 
observed (see table 3.13 of Biretta et al 2000). This is 
the reason why, for example, as shown in Table 8, the ex- 
tinction coefficient TZ in the F170W filter is about 4 times 
smaller for a Te// = 3,500 K star than for a 40,000 K 
one! One can easily assess the influence of the red leak by 
recomputing the same coefficients after cutting it out from 
the response curve of the filter, in this case the F170W. 
What one finds, then, is that, while the extinction coef- 
ficient for Te// ~ 40,000 K is virtually unchanged, the 
one for T,.// = 3. 500 K increases by a factor of 3, thus 
reducing the difference between them to about 30%. The 
residual difference in H, between different temperatures, 
the lower one yielding to lower extinction coefficients, is 
caused by the fact that filters have a finite, non negligible 
bandwidth and, thus, the effective extinction in a band 
depends on the shape of the spectrum which is a function 
of temperature. The same explanation applies for the fact 
(also illustrated in Table 8) that extinction coefficients be- 
come smaller as the optical depth of the dust increases 
because dust extinction makes a spectral energy distribu- 
tion become redder. 

Tables with similar quantities computed for all of 
WFPC2 filters and for all model atmospheres in the Bessel 
et al (1998) grid are available and may be provided on re- 
quest. 
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See fig01.jpg 

Fig. 1. — The field centered on SN1987A (about 130" radius) as observed in the combination of the B, V, and I broad bands plus the [OIII] 
and Ha narrow band images. 



See fig02.jpg 

Fig. 2. — Color-Magnitude Diagrams for four combination of filters. The grey dots are stars with average errors, as defined in equation (1), 
i5b > 0.1, whereas the black dots are the 6,695 stars with Ss < 0.1. The reddening vector corresponding to E^B — V) = 0.2 is show in each 
panel 



Fig. 3.— 



See fig03.jpg 

Photometric error as a function of magnitude for the 6 broad band filters around SN1987A. The 0.1 mag error level is indicated. 
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Fig. 4. — Four examples of the spectral resolution of the "Wide Band spectroscopy". The magnitudes and errors of these four stars are 
given in Table 3. The position of the Balmer jump at 3646 A is indicated with a dashed line. 
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Fig. 5. — Comparison of spectra for different combinations of Tg/y and B(B — V): 
T^ff = 7,600 and S(B - V) = 0.55 (long-dash line), T^yj = 12,000 and E(B - V) = 

F814W flux and errorbars corresponding to 0.1 mag are shown. As it can be seen, at this resolution the spectra are almost indistinguishable 
for wavelengths longer than 3,300 A and ultraviolet imaging is needed to tell them apart. 



Te// = 5, 000 and E(B - V) = (short-dash Une), 
0.86 (full line). The spectra are normalized to the 
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Fig. 6.— Reddening-free color QuBi as a function of (F336W - F814W) from the theoretical models by Bessel et al (1998) for Z= 0.3-Z© 
and log(g) = 4.5. The horizontal arrow indicates the reddening vector for £(B — V) = 0.2 mag and a few important temperatures along the 
sequence are also marked. Stars in the four regions A — » £) have a different number of solutions (see text). The dashed line highlights the 
location in this diagram of the T Tauri stars: their colors axe not well represented by normal photospheres. 



Accurate stellar population studies 




7. — Location of class I (filled circles), II (filled triangles) and III (small dots) stars in the QuBi vs. F336W — F814W plane. 
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Fig. 8. — Same as Figure 4, but with the best fitting model shown as a full line. The parameters derived from the fit are listed in Table 4. 



See fig09.jpg 



Fig. 9. — HR diagram for the stars in the field around SN1987A. The black dots are the 9,474 stars with <51og(Te//) < 0.05. The Zero 
Age Main Sequence for Z= 0.3 • Zq from the models by Brocato &c Castellani (1993) and Cassisi et al (1994) for masses below 25 Mq and by 
Schaerer et al (1993) above 25 Mq is shown as a full line. 
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See figlOa.jpg 

Fig. 10. — Results of the "blind" fit, for the stars around SN1987A: the fit described in section 4.2 is performed directly without the 
selection of section 4.1. Left panel: HR diagram The grey triangles highlight the stars for which the fit returns too high a value of T^eff ^'id 
E(B — V). These are stars with multiple solutions for which the "blind" fit picks the wrong one (see text). The Zero Age Main Sequence for 
Z= 0.3 • Zq from the models by Brocato & Castellani (1993) and Cassisi et al (1994) for masses below 25 M© and by Schaerer et al (1993) 
above 25 Mq is shown as a full line. Right panel: E(B — V) distribution resulting from the "blind" fit. The high reddening peak (hatched 
histogram) is entirely made of stars plotted in grey in the left panel, thus confirming that their location in the HR diagram is due to an 
erroneous reddening determination. 
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Fig. 11. — Panel (a): reddening histogram from the stars in the field of Supernova 1987A. Only stars for which E(B — V) was measured 
individually (the 2,510 class I and II stars, see the text) are used here. Panel (b): spatial distribution of stars with high {E(B — V) > 0.3, 
filled circles) and low (B(B — V) < 0.1, open circles) reddening. 



20 



Romaniello et al 



See figl2.jpg 

Fig. 12. — Before and after the reddening correction. The observed F439W vs. F439W — F814W diagram for all the stars in our WFPC2 
frames with 5^ < 0.1 is shown in panel (a), while the dereddened one in panel (b). Black dots represent the stars for which E(B — V) was 
determined individually (class I and II, see text), grey dots stars dereddened with the mean value of E(B — V) from the closest class I and II 
neighbors. The density contours, spaced by factors of 2, are are superposed in red. The visual impression that the Main Sequence in panel (b) 
is narrower than the one in panel (a) is confirmed by the histograms shown in panel (c). Here, the blue line represents the color distribution 
of the dereddened stars in the 21 < m(F439W)o < 22 magnitude interval and the red one the distribution of the observed colors for the same 
stars. For the sake of clarity, the latter was shifted by 0.5 mag in color to have the same mode of the dereddened one. There are 218 and 192 
stars, respectively, at the peak of the dereddened and observed distribution, while both histograms contain a total of 1923 stars. 
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Fig. 13. — Width of the Main Sequence as a function of magnitude before (red line) and after (blue line) reddening correction. The 
magnitude interval along the Main Sequence is indicated in each panel. The faintest magnitude bin heavily suffers from incompleteness. Also, 
the raw histograms were shifted in color in order to have the same mode as the corrected ones. As it can be seen, the corrected histograms 
are systematically narrower than the uncorrected ones. 
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See figl4.jpg 

Fig. 14. — Same as Figure 12, but the contours in panels (a) and (b) and the histograms in panel (c) refer only to the stars with individual 
reddening correction (class I and II, see text). The peak values of the histograms in panel (c) are 174 and 141 for the corrected (blue line) 
and observed one (red line), respectively. 
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Fig. 15. — Zcro-rcddcning QuBi vs. F336W — F814W relation for three values of surface gravity: log(g) = 5 (dashed line), log(g) = 4.5 
(solid line) and log(g) = 4 (dotted line). The model atmospheres are those of Bessel et al (1998) for Z= 0.3 • Zq. The points for T^f j=50,000, 
15,000, 9,250, 6,000 and 3,500 K are marked on each curve. 
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Fig. 16. — Zero-rcddcniiig Qvbi vs. F336W - F814W relation for three values of metallicity: Z=Zq (dashed line), Z= 0.3 • Zq (sohd hne) 
and Z=Zq/10 (dotted line). The model atmospheres are those of Bessel et al (1998) for log(g) = 4.5. The points for Tejj=50,000, 15,000, 
9,250, 6,000 and 3,500 K are marked on each curve. 
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Fig. 17.— The Qbvi reddening-free color as a function of F439W - F814W (panel (a)) and Te// (panel (b)) for the Bessel et al (1998) 
model atmospheres for Z= 0.3-Zq and log{g) = 4.5. The case of black body spectra is shown in both panels as a dashed line. 
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Fig. 18. — Contour plots in the Te//-i7(B — V) space for four stars with unambiguous solution. Left panel: only four bands (F439W, 
F555W, F675W and F814W) are used in the fit. Right panel: all six bands are used. In both panels the contours are for Sx^ = 2 from the 



minimum. 
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Fig. 20.— contour map for artificial stars with E(B - V) = 0.2 and Te// = 20,000 K (panel (a)), 25,000 K (b), 30,000 K (c) and 
40,000 K (d). The contours correspond to the Icr confidence level for a 5-band fit. In all cases the same 4 optical filters arc used (F439W 
(B), F555W (V), F675W (R) and F814W (I)) in combination with difl'erent UV passbands; F170W (dot-dashed line), F255W (dashed line) 
or F336W (full line). The photometric error is assumed to be the same in all 3 UV filters. The thin straight lines mark the input parameters 
of the model stars. 
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See fig21.jpg 

Fig. 21.— HR diagram for the 13,098 stars detected in the control field for SN1987A. Black dots arc the 4,912 stars with <51og{Te//) < 0.05. 
The open squares indicate stars for which the F300W magnitude is ill determined because of saturation and, hence, the fit was performed 
excluding this filter. The location of Sk -69 211 according to the photometry by Fitzpatrick (1998) is shown with a star symbol. The Zero 
Age Main Sequence for Z= 0.3 • Zq from the models by Brocato & Castellani (1993) and Cassisi et al (1994) for masses below 25 Mq and by 
Schaerer et al (1993) above 25 Mq is shown as a full line. 



See fig22.jpg 

Fig. 22. — Ratio of output to input temperature as a function of the latter for 10,000 model stars distributed evenly in log(Tin). The fit 
was performed as described in section 4 using 5 bands: F336W, F439W, F555W, F675W and F814W. The input stars in panel (a) simulate a 
random photometric error of 0.1 mag in all of the bands, those in panel (b) an error of 0.05 mag and those in panel (c) an error of 0.02 mg. 
Black dots are class I and II stars for which the reddening was determined individually, while grey ones are class III stars dereddened with 
the average value of their class I and II neighbors. 



See fig23.jpg 

Fig. 23. — Same as Figure 22, but the fit was performed with different combinations of bands: 4 in panel (a) (F336W, F439W, F555W and 
F814W), 5 in panel (b) (those of panel (a) plus F675W) and 6 in panel (c) (those of panel (b) plus F255W). The stars simulate a random 
photometric error of 0.05 mag. 
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Fig. 24. — Distribution of the ratio of output to input temperature in the case of systematic errors for four representative temperatures: 
5,000 K in panel (a), 10,000 K in panel (b), 20,000 K in panel (c) and 40,000 K in panel (d). In each panel the dotted line is for the case of a 
systematic error of 0.02 mag, the full one for 0.05 mag and the dashed one for 0.1 mag. 
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Table 1 

Log of the observations centered on Supernova 1987A. 



Filter Name 



Exposure Time (seconds) 



Comments 



September 1994^^ February 1996'' July 1997^= 



F255W 
F336W 
F439W 
F555W 
F675W 
F814W 



2x900 
2x600 
2x400 
2x300 
2x300 
2x300 



1100+1400 

2x600 
350+600 
2x300 
2x300 
2x300 



2x1300 
2x800 
2x400 
2x300 
2x300 
2x400 



UV Filter 
U Filter 
B Filter 
V Filter 
R Filter 
I Filter 



"September 24, 1994, proposal number 5753. 
''February 6, 1996, proposal number 6020. 
°July 10, 1997, proposal number 6437. 



Table 2 

Log of the observations of the control field for SN1987A. They were taken on January, 1"* 1997 and the 

PROPOSAL NUMBER IS 6437. 

Filter Name Exposure Time (s) Comments 

F300W 600 + 1200 Wide U filter 

F450W 50 + 200 Wide B filter 

F675W 50 + 200 R filter 

F814W 140 + 400 I filter 



Table 3 

Observed magnitudes and statistical errors for the stars of Figure 4. 



Panel F255W 


F336W 


F439W 


F555W 


F675W 


F814W 


(a) 15.74 + 0.01 

(b) 19.34 ± 0.04 

(c) 21.50 + 0.20 

(d) 22.10 + 0.30 


16.13 + 0.01 
19.24 + 0.01 
20.99 + 0.03 
21.32 + 0.04 


17.37 + 0.01 
19.83 + 0.01 
20.91 + 0.02 
21.23 + 0.02 


17.23 + 0.01 
19.69 + 0.01 
20.71 + 0.01 
20.85 + 0.01 


17.44 + 0.01 

19.45 + 0.01 

20.45 + 0.01 

20.46 + 0.01 


17.48 + 0.01 
19.28 + 0.01 
20.35 + 0.01 
20.26 + 0.01 



Table 4 

Best-fit parameters for the spectra of Figure 8. The star in panel (a) is saturated and this causes the large 
quoted errors on the quantities derived for it. The corresponding magnitudes are listed in Table 3. 

Panel Te//(K) E{B~\) log(L/L0) R/Rq 

(a) 25, 350 + 5, 700 0.202 + 0.080 3.662 + 0.264 3.444 + 0.090 

(b) 18,160 + 2,100 0.431 + 0.051 2.695 + 0.166 2.265 + 0.050 

(c) 8,000 + 350 0.117 + 0.050 1.176 + 0.079 2.070 + 0.042 

(d) 7,480 + 300 0.209 + 0.047 1.255 + 0.067 2.548 + 0.052 
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Table 5 

Mean ratio and rms scatter of the output to input values of selected quantities for 10,000 artificial stars 

(see text). 



Input 


Fitted 


Te 


ff 




-V) 


L 




error 
(rms) 


bands 


Tin < 9, 000 


Tin > 9, 000 


Tin < 9, 000 


Tin > 9, 000 


Tin < 9, 000 


Tin > 9, 000 


0.1 
0.05 
0.02 
0.05 
0.05 


5 
5 
5 
4 
6 


0.99 ±0.05 
0.99 ±0.03 
0.99 ±0.04 
0.99 ± 0.04 
1.00 ± 0.03 


0.97 ±0.18 
0.98 ±0.12 
0.99 ±0.06 
0.98 ±0.12 
0.98 ±0.11 


0.01 ±0.02 
0.01 ±0.02 
0.01 ±0.03 
0.01 ± 0.02 
0.01 ±0.02 


-0.02 ±0.06 
-0.01 ±0.04 
-0.00 ±0.02 
-0.01 ±0.04 
-0.01 ±0.04 


0.98 ±0.17 
0.97 ±0.10 
0.97 ±0.11 
0.98 ±0.11 
0.99 ±0.10 


1.01 ±0.58 
0.99 ±0.37 
0.97 ±0.18 
0.99 ±0.38 
0.97 ±0.35 
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Table 6 

Magnitudes in selected WFPC2 broad band filters from the Bessel et al (1998) models for Z=Zq = 0.02 

REFERRED TO A 1 Rq STAR AT A DISTANCE OF 10 PC (SEE EQUATION (20)). 
Te// F170W F255W F300W F336W F439W F450W F555W F606W F675W F814W 









3500 


11.930 


16.528 


4000 


10.692 


14.882 


4500 


9.900 


12.670 


5000 


9.250 


9.935 


5500 


8.665 


8.278 


6000 


8.085 


7.105 


6500 


7.381 


6.183 


7000 


6.501 


5.442 


7500 


5.647 


4.829 


8000 


4.836 


4.274 


8500 


4.077 


3.745 


9000 


3.415 


3.278 


9500 


2.886 


2.869 


10000 


2.456 


2.503 


10500 


2.089 


2.178 


11000 


1.768 


1.886 


11500 


1.480 


1.622 


12000 


1.220 


1.381 


12500 


0.981 


1.161 


13000 


0.762 


0.959 


15000 


0.034 


0.290 


17000 


-0.537 


-0.242 


19000 


-1.016 


-0.695 


21000 


-1.441 


-1.091 


23000 


-1.832 


-1.435 


25000 


-2.197 


-1.737 


27000 


-2.492 


-1.918 








3500 


11.574 


15.806 


4000 


10.653 


14.317 


4500 


9.904 


12.589 


5000 


9.250 


9.912 


5500 


8.654 


8.115 


6000 


8.047 


6.844 


6500 


7.301 


5.855 


7000 


6.337 


5.065 


7500 


5.459 


4.433 


8000 


4.703 


3.913 


8500 


3.980 


3.472 


9000 


3.330 


3.082 


9500 


2.801 


2.731 


10000 


2.364 


2.417 


10500 


1.999 


2.128 


11000 


1.682 


1.863 


11500 


1.401 


1.620 


12000 


1.146 


1.396 


12500 


0.914 


1.190 


13000 


0.700 


0.999 


15000 


-0.016 


0.361 


17000 


-0.577 


-0.137 


19000 


-1.040 


-0.550 


21000 


-1.440 


-0.909 


23000 


-1.800 


-1.231 


25000 


-2.132 


-1.525 


27000 


-2.446 


-1.792 


29000 


-2.744 


-2.040 


31000 


-3.015 


-2.271 


33000 


-3.238 


-2.472 


35000 


-3.412 


-2.632 


37000 


-3.551 


-2.755 


39000 


-3.666 


-2.851 


41000 


-3.768 


-2.933 


43000 


-3.862 


-3.008 


45000 


-3.949 


-3.077 
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-4.030 


-3.142 
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1.913 
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1.581 
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1.724 


1.345 


1.393 


1.461 


1.557 


1.165 


1.220 


1.295 


1.403 


0.995 


1.054 


1.137 


1.254 


0.831 


0.895 


0.984 


1.109 


0.676 


0.743 


0.839 


0.971 


0.527 


0.598 


0.699 


0.838 


0.378 


0.452 


0.558 


0.704 


0.237 


0.313 


0.423 


0.572 


0.117 


0.193 


0.304 


0.456 


0.021 


0.097 


0.207 


0.360 


-0.053 


0.023 


0.132 


0.286 


-0.115 


-0.039 


0.071 


0.226 


-0.173 


-0.096 


0.015 


0.171 


-0.228 


-0.150 


-0.039 


0.118 


-0.281 


-0.202 


-0.090 


0.067 
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Table 6 — Continued 





F170W 


F255W 


F300W 


F336W 


F439W 


F450W 


F555W 


F606W 


F675W 


F814W 


49000 


-4.106 


-3.202 


-2.727 


-2.323 


-0.595 


-0.570 


-0.331 


-0.251 


-0.139 


0.019 



Table 7 

Magnitudes in selected WFPC2 broad band filters from the Bessel et al (1998) models for Z=O.3Z0 = 6 • 10~^ 

REFERRED TO A 1 R© STAR AT A DISTANCE OF 10 PC (SEE EQUATION (20)). 



F170W F255W F300W F336W F439W F450W F555W F606W F675W F814W 



'-eff 



log(g) = 3 




3500 


11.592 


4000 


10.663 


4500 


9.877 


5000 


9.212 


5500 


8.586 


6000 


7.880 


6500 


7.004 


7000 


6.157 


7500 


5.395 


8000 


4.637 


8500 


3.927 


9000 


3.318 


9500 


2.825 


10000 


2.414 


11000 


1.748 


12000 


1.213 


13000 


0.764 


15000 


0.039 


17000 


-0.538 


19000 


-1.025 


21000 


-1.452 


23000 


-1.839 


25000 


-2.199 


log(g) = 4.5 




3500 


11.438 


4000 


10.589 


4500 


9.883 


5000 


9.214 


5500 


8.570 


6000 


7.827 


6500 


6.882 


7000 


5.920 


7500 


5.158 


8000 


4.490 


8500 


3.827 


9000 


3.207 


9500 


2.724 


10000 


2.321 


11000 


1.673 


12000 


1.154 


13000 


0.716 


15000 


0.005 


17000 


-0.560 


19000 


-1.031 


21000 


-1.439 


23000 


-1.804 


25000 


-2.138 


27000 


-2.450 


29000 


-2.744 


31000 


-3.015 


33000 


-3.242 


35000 


-3.418 


40000 


-3.716 



16.214 


12.581 


11.833 


14.631 


11.027 


10.147 


11.619 


9.048 


8.245 


9.214 


7.539 


6.909 


7.704 


6.430 


5.910 


6.595 


5.617 


5.189 


5.761 


5.035 


4.694 


5.119 


4.591 


4.329 


4.590 


4.211 


4.017 


4.114 


3.848 


3.711 


3.646 


3.478 


3.394 


3.214 


3.123 


3.085 


2.819 


2.788 


2.788 


2.465 


2.479 


2.508 


1.872 


1.945 


2.013 


1.392 


1.510 


1.611 


0.988 


1.150 


1.283 


0.339 


0.574 


0.765 


-0.183 


0.107 


0.342 


-0.627 


-0.296 


-0.027 


-1.016 


-0.655 


-0.360 


-1.359 


-0.974 


-0.659 


-1.670 


-1.266 


-0.934 


15.667 


12.222 


11.409 


14.012 


10.782 


9.907 


11.779 


9.151 


8.312 


9.190 


7.549 


6.889 


7.566 


6.353 


5.824 


6.369 


5.446 


5.016 


5.421 


4.729 


4.378 


4.741 


4.271 


4.015 


4.190 


3.881 


3.698 


3.738 


3.548 


3.426 


3.355 


3.251 


3.175 


3.018 


2.975 


2.935 


2.692 


2.702 


2.698 


2.388 


2.441 


2.468 


1.853 


1.966 


2.042 


1.406 


1.560 


1.671 


1.024 


1.211 


1.349 


0.404 


0.646 


0.833 


-0.087 


0.203 


0.432 


-0.497 


-0.168 


0.094 


-0.854 


-0.494 


-0.204 


-1.172 


-0.786 


-0.473 


-1.461 


-1.054 


-0.721 


-1.726 


-1.300 


-0.951 


-1.975 


-1.533 


-1.167 


-2.212 


-1.758 


-1.379 


-2.425 


-1.964 


-1.578 


-2.597 


-2.134 


-1.744 


-2.868 


-2.401 


-2.009 



10.214 


9.855 


8.598 


8.732 


8.331 


7.282 


7.504 


7.120 


6.323 


6.576 


6.261 


5.621 


5.811 


5.576 


5.074 


5.168 


5.001 


4.620 


4.622 


4.504 


4.218 


4.147 


4.064 


3.856 


3.722 


3.670 


3.533 


3.354 


3.331 


3.263 


3.078 


3.068 


3.045 


2.870 


2.865 


2.875 


2.704 


2.702 


2.733 


2.565 


2.564 


2.609 


2.331 


2.332 


2.396 


2.135 


2.137 


2.215 


1.964 


1.967 


2.058 


1.671 


1.677 


1.791 


1.416 


1.425 


1.558 


1.185 


1.196 


1.344 


0.972 


0.983 


1.144 


0.774 


0.786 


0.956 


0.583 


0.594 


0.775 


10.067 


9.678 


8.397 


8.693 


8.320 


7.249 


7.496 


7.163 


6.342 


6.569 


6.267 


5.625 


5.839 


5.595 


5.079 


5.231 


5.048 


4.641 


4.743 


4.602 


4.275 


4.279 


4.179 


3.930 


3.885 


3.810 


3.615 


3.531 


3.478 


3.329 


3.217 


3.186 


3.084 


2.966 


2.953 


2.893 


2.771 


2.767 


2.738 


2.610 


2.612 


2.608 


2.356 


2.363 


2.395 


2.154 


2.163 


2.218 


1.981 


1.991 


2.062 


1.690 


1.703 


1.800 


1.446 


1.462 


1.580 


1.231 


1.249 


1.385 


1.034 


1.054 


1.206 


0.851 


0.871 


1.037 


0.677 


0.698 


0.876 


0.511 


0.533 


0.722 


0.349 


0.372 


0.571 


0.186 


0.210 


0.420 


0.031 


0.055 


0.273 


-0.101 


-0.077 


0.147 


-0.319 


-0.296 


-0.064 



8.116 


7.605 


6.648 


6.832 


6.354 


5.809 


5.973 


5.585 


5.166 


5.338 


5.016 


4.675 


4.841 


4.570 


4.292 


4.431 


4.206 


3.980 


4.071 


3.891 


3.713 


3.747 


3.613 


3.480 


3.463 


3.373 


3.281 


3.225 


3.174 


3.114 


3.028 


3.004 


2.969 


2.869 


2.861 


2.844 


2.734 


2.736 


2.731 


2.614 


2.622 


2.628 


2.406 


2.422 


2.444 


2.231 


2.253 


2.288 


2.079 


2.108 


2.154 


1.820 


1.862 


1.927 


1.595 


1.648 


1.729 


1.388 


1.451 


1.545 


1.193 


1.262 


1.367 


1.009 


1.084 


1.196 


0.831 


0.912 


1.032 


7.860 


7.299 


6.488 


6.775 


6.273 


5.722 


5.979 


5.575 


5.158 


5.340 


5.013 


4.675 


4.844 


4.569 


4.290 


4.446 


4.214 


3.982 


4.113 


3.916 


3.722 


3.802 


3.646 


3.495 


3.519 


3.400 


3.288 


3.263 


3.181 


3.105 


3.046 


2.998 


2.954 


2.876 


2.854 


2.833 


2.734 


2.729 


2.726 


2.613 


2.619 


2.629 


2.410 


2.430 


2.458 


2.239 


2.268 


2.308 


2.088 


2.124 


2.175 


1.834 


1.880 


1.949 


1.621 


1.678 


1.761 


1.433 


1.499 


1..594 


1.259 


1.333 


1.438 


1.095 


1.177 


1.291 


0.939 


1.027 


1.149 


0.788 


0.882 


1.012 


0.642 


0.741 


0.878 


0.493 


0.598 


0.741 


0.349 


0.457 


0.605 


0.223 


0.331 


0.481 


0.011 


0.116 


0.267 
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Table 8 

Extinction coefficients 72. in selected WFPC2 broad band filters from the Bessel et al (1998) model 

ATMOSPHERES AND THE GALACTIC REDDENING LAW OF SAVAGE & MATHIS (1979) AS A FUNCTION OF EFFECTIVE TEMPERATURE 

AND E{B — V)(SEE EQUATION (20)). 



£;(B-V) Te//=3500 K 4000 5250 6500 8750 9250 10000 25000 40000 



F170W 








0.0 


2.017 


2.205 


3.060 


0.2 


1.970 


2.143 


2.851 


0.4 


1.927 


2.088 


2.689 


0.6 


1.889 


2.039 


2.559 


0.8 


1.854 


1.994 


2.452 


1.0 


1.822 


1.954 


2.361 


1.2 


1.793 


1.916 


2.284 


1.4 


1.766 


1.882 


2.216 


1.5 


1.753 


1.866 


2.185 


F255W 








0.0 


4.699 


5.262 


6.256 


0.2 


4.437 


5.095 


6.214 


0.4 


4.147 


4.892 


6.171 


0.6 


3.848 


4.652 


6.125 


0.8 


3.562 


4.385 


6.072 


1.0 


3.304 


4.108 


6.005 


1.2 


3.081 


3.840 


5.914 


1.4 


2.892 


3.595 


5.787 


1.5 


2.809 


3.483 


5.707 


F300W 








0.0 


3.789 


4.339 


5.173 


0.2 


3.577 


4.185 


5.131 


0.4 


3.360 


4.004 


5.079 


0.6 


3.152 


3.804 


5.014 


0.8 


2.965 


3.598 


4.931 


1.0 


2.802 


3.399 


4.822 


1.2 


2.664 


3.217 


4.685 


1.4 


2.547 


3.055 


4.520 


1.5 


2.496 


2.982 


4.429 


F336W 








0.0 


4.045 


4.412 


4.911 


0.2 


3.888 


4.296 


4.888 


0.4 


3.721 


4.160 


4.857 


0.6 


3.554 


4.008 


4.816 


0.8 


3.397 


3.848 


4.761 


1.0 


3.254 


3.689 


4.690 


1.2 


3 1 30 


3 53Q 


4 5QQ 


1.4 


3.023 


3.403 


4.489 


1.5 


2.976 


3.341 


4.428 


F439W 








0.0 


4.088 


4.097 


4.132 


0.2 


4.085 


4.094 


4.129 


0.4 


4.083 


4.092 


4.127 


0.6 


4.080 


4.089 


4.124 


0.8 


4.076 


4.086 


4.121 


1.0 


4.073 


4.083 


4.118 


1.2 


4.069 


4.079 


4.116 


1.4 


4.065 


4.076 


4.113 


1.5 


4.062 


4.074 


4.111 


F450W 








0.0 


3.787 


3.798 


3.835 



6, 


.202 


7, 


.901 


7, 


.956 


8, 


.009 


8, 


.110 


8, 


.111 


5, 


.484 


7, 


.717 


7, 


.818 


7, 


.908 


8, 


.080 


8, 


.086 


4, 


.812 


7, 


.427 


7, 


.598 


7, 


.746 


8, 


.045 


8, 


.058 


4, 


.265 


6 


.993 


7, 


.248 


7, 


.480 


7, 


.998 


8, 


.023 


3, 


.849 


6, 


.438 


6, 


.757 


7, 


.073 


7, 


.926 


7, 


.973 


3, 


.534 


5, 


.856 


6, 


.192 


6, 


.550 


7, 


.799 


7, 


.892 


3 


.289 


5, 


.331 


5, 


.649 


6, 


.003 


7, 


.578 


7, 


.745 


3, 


.095 


4, 


.893 


5, 


.181 


5, 


.508 


7, 


.232 


7, 


.491 


3, 


.012 


4, 


.704 


4, 


.977 


5, 


.288 


7, 


.021 


7, 


.319 


6, 


.548 


6, 


.772 


6, 


.791 


6, 


.812 


6, 


.934 


6, 


.958 


6, 


.494 


6 


.713 


6, 


.731 


6, 


.751 


6 


.871 


6, 


,895 


6, 


.443 


6, 


.657 


6, 


.675 


6, 


.696 


6, 


.813 


6, 


.836 


6, 


,396 


6. 


.606 


6, 


.623 


6, 


.644 


6, 


.759 


6, 


,781 


6, 


.350 


6 


.557 


6, 


.575 


6, 


.595 


6 


.709 


6, 


,731 


6, 


.303 


6 


.510 


6, 


.528 


6, 


.549 


6 


.662 


6, 


,684 


6, 


.252 


6, 


.463 


6, 


.482 


6, 


.503 


6, 


.619 


6, 


,639 


6, 


.192 


6, 


.414 


6, 


.434 


6, 


.457 


6, 


.577 


6. 


,598 


6, 


.156 


6, 


.388 


6, 


.408 


6, 


.433 


6, 


.557 


6, 


,578 


5, 


.391 


5, 


.603 


5, 


.628 


5, 


.660 


5, 


.851 


5, 


,883 


5, 


.353 


5, 


.558 


5, 


.582 


5, 


.614 


5, 


.804 


5, 


,836 


5, 


.314 


5, 


.513 


5, 


.537 


5, 


.569 


5, 


.759 


5, 


,791 


5, 


.271 


5, 


.467 


5, 


.492 


5, 


.524 


5, 


.716 


5, 


,747 


5, 


.221 


5, 


.419 


5, 


.444 


5, 


.477 


5, 


.674 


5, 


,705 


5, 


.160 


5, 


.365 


5, 


.391 


5, 


.427 


5, 


.632 


5, 


,665 


5, 


.083 


5, 


.301 


5, 


.330 


5, 


.370 


5, 


.591 


5, 


,624 


4, 


.982 


5, 


.222 


5, 


.255 


5, 


.301 


5, 


.549 


5, 


,584 


4, 


.921 


5, 


.175 


5, 


.211 


5, 


.260 


5, 


.526 


5, 


,563 


4, 


.974 


4, 


.998 


5, 


.002 


5, 


.007 


5, 


.037 


5. 


,041 


4, 


.962 


4, 


.990 


4, 


.994 


5, 


.000 


5. 


.033 


5, 


,037 


4, 


.948 


4, 


.979 


4, 


.984 


4, 


.991 


5, 


.028 


5, 


,033 


4, 


.929 


4, 


.965 


4, 


.971 


4, 


.979 


5, 


.022 


5, 


,028 


4, 


.903 


4, 


.947 


4, 


.954 


4, 


.964 


5, 


.016 


5, 


,022 


4, 


.868 


4, 


.922 


4, 


.931 


4, 


.944 


5, 


.009 


5, 


,016 


4, 


.821 


4, 


.888 


4, 


.900 


4, 


.917 


5, 


.000 


5, 


,009 


4, 


.759 


4, 


.842 


4, 


.858 


4, 


.880 


4, 


.988 


4, 


,999 


4, 


.721 


4, 


.814 


4, 


.831 


4, 


.857 


4, 


.980 


4, 


,994 


4, 


.152 


4, 


.161 


4, 


.163 


4, 


.166 


4, 


.178 


4, 


,181 


4, 


.149 


4, 


.159 


4, 


.161 


4, 


.164 


4, 


.176 


4, 


,179 


4, 


.147 


4, 


.156 


4, 


.158 


4, 


.161 


4, 


.173 


4, 


,177 


4, 


.144 


4, 


.154 


4, 


.156 


4, 


.159 


4, 


.171 


4, 


,174 


4, 


.142 


4, 


.152 


4, 


.154 


4, 


.156 


4, 


.168 


4, 


,172 


4, 


.139 


4, 


.149 


4, 


.151 


4, 


.154 


4, 


.166 


4, 


,169 


4, 


.136 


4, 


.147 


4, 


.149 


4, 


.151 


4, 


.163 


4, 


,167 


4, 


.134 


4, 


.144 


4, 


.146 


4, 


.149 


4, 


.161 


4, 


,164 


4, 


.132 


4, 


.143 


4, 


.145 


4, 


.148 


4, 


.160 


4, 


,163 


3, 


.890 


3 


.929 


3, 


.935 


3, 


.944 


3. 


.987 


4, 


,000 
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Table 8 — Continued 



E{B - V) 


Teff=3500 K 


4000 


5250 


6500 


8750 


9250 


10000 


25000 


40000 


0.2 


3.779 


3.790 


3.825 


3.880 


3.918 


3.925 


3.933 


3.976 


3.989 


0.4 


3.772 


3.782 


3.816 


3.869 


3.908 


3.914 


3.922 


3.965 


3.977 


0.6 


3.764 


3.774 


3.806 


3.859 


3.897 


3.903 


3.912 


3.953 


3.966 


0.8 


3.757 


3.767 


3.797 


3.848 


3.887 


3.893 


3.901 


3.942 


3.954 


1.0 


3.749 


3.759 


3.789 


3.838 


3.876 


3.882 


3.890 


3.930 


3.942 


1.2 


3.742 


3.752 


3.780 


3.828 


3.866 


3.872 


3.880 


3.919 


3.931 


1.4 


3.735 


3.745 


3.771 


3.818 


3.856 


3.862 


3.869 


3.908 


3.919 


1.5 


3.731 


3.741 


3.767 


3.814 


3.851 


3.857 


3.864 


3.902 


3.914 


F555W 




















0.0 


3.031 


3.066 


3.153 


3.189 


3.243 


3.249 


3.255 


3.281 


3.288 


0.2 


3.019 


3.053 


3.139 


3.175 


3.228 


3.235 


3.241 


3.267 


3.274 


0.4 


3.008 


3.041 


3.126 


3.162 


3.214 


3.220 


3.227 


3.252 


3.259 


0.6 


2.997 


3.030 


3.113 


3.148 


3.200 


3.206 


3.212 


3.238 


3.245 


0.8 


2.987 


3.018 


3.100 


3.135 


3.186 


3.193 


3.199 


3.224 


3.230 


1.0 


2.977 


3.007 


3.088 


3.122 


3.173 


3.179 


3.185 


3.210 


3.216 


1.2 


2.967 


2.997 


3.076 


3.110 


3.159 


3.166 


3.172 


3.196 


3.203 


1.4 


2.958 


2.987 


3.064 


3.097 


3.146 


3.153 


3.158 


3.183 


3.189 


1.5 


2.953 


2.982 


3.058 


3.091 


3.140 


3.146 


3.152 


3.176 


3.182 


F606W 




















0.0 


2.712 


2.729 


2.805 


2.841 


2.893 


2.900 


2.906 


2.932 


2.938 


0.2 


2.703 


2.719 


2.793 


2.828 


2.879 


2.886 


2.892 


2.918 


2.924 


0.4 


2.694 


2.709 


2.781 


2.815 


2.866 


2.873 


2.879 


2.904 


2.910 


0.6 


2.685 


2.700 


2.770 


2.803 


2.853 


2.860 


2.866 


2.891 


2.897 


0.8 


2.677 


2.691 


2.759 


2.791 


2.841 


2.847 


2.853 


2.877 


2.883 


1.0 


2.669 


2.682 


2.748 


2.780 


2.828 


2.834 


2.840 


2.864 


2.870 


1.2 


2.661 


2.674 


2.737 


2.768 


2.816 


2.822 


2.828 


2.852 


2.857 


1.4 


2.654 


2.665 


2.727 


2.758 


2.805 


2.810 


2.816 


2.839 


2.844 


1.5 


2.650 


2.661 


2.722 


2.752 


2.799 


2.805 


2.810 


2.833 


2.838 


F675W 




















0.0 


2.430 


2.432 


2.444 


2.450 


2.460 


2.461 


2.462 


2.468 


2.470 


0.2 


2.427 


2.429 


2.441 


2.447 


2.457 


2.458 


2.460 


2.465 


2.467 


0.4 


2.425 


2.426 


2.439 


2.445 


2.454 


2.456 


2.457 


2.463 


2.464 


0.6 


2.422 


2.424 


2.436 


2.442 


2.452 


2.453 


2.454 


2.460 


2.461 


0.8 


2.419 


2.421 


2.433 


2.439 


2.449 


2.450 


2.451 


2.457 


2.459 


1.0 


2.416 


2.418 


2.430 


2.436 


2.446 


2.447 


2.449 


2.455 


2.456 


1.2 


2.413 


2.416 


2.428 


2.434 


2.443 


2.445 


2.446 


2.452 


2.453 


1 zL 


Z. 41:1 1 




Z.41Z0 


Z.4101 


9 zLzLI 

Z. 41:41:1 


9 AAO 

Z.41:41:Z 


9 AA'i 
Z. 41410 


9 AAQ 




1.5 


2.409 


2.412 


2.424 


2.430 


2.439 


2.441 


2.442 


2.448 


2.449 


F814W 




















0.0 


1.786 


1.829 


1.860 


1.875 


1.892 


1.894 


1.896 


1.912 


1.916 


0.2 


1.777 


1.820 


1.851 


1.867 


1.884 


1.886 


1.888 


1.904 


1.908 


0.4 


1.768 


1.811 


1.842 


1.858 


1.876 


1.877 


1.879 


1.896 


1.900 


0.6 


1.759 


1.801 


1.833 


1.849 


1.867 


1.868 


1.870 


1.888 


1.892 


0.8 


1.750 


1.792 


1.824 


1.840 


1.858 


1.859 


1.861 


1.879 


1.884 


1.0 


1.740 


1.782 


1.815 


1.831 


1.849 


1.850 


1.852 


1.870 


1.875 


1.2 


1.731 


1.773 


1.805 


1.822 


1.839 


1.841 


1.843 


1.861 


1.866 


1.4 


1.721 


1.763 


1.795 


1.812 


1.830 


1.832 


1.833 


1.852 


1.857 


1.5 


1.717 


1.758 


1.791 


1.807 


1.825 


1.827 


1.829 


1.848 


1.853 
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